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Abstract. Space-charge dominated streamer discharges can emerge in free space 
from single electrons. We reinvestigate the Raether-Meek criterion and show that 
streamer emergence depends not only on ionization and attachment rates and gap 
length, but also on electron diffusion. Motivated by simulation results, we derive 
an explicit quantitative criterion for the avalanche-to-streamer transition both for 
pure non-attaching gases and for air, under the assumption that the avalanche 
emerges from a single free electron and evolves in a homogenous field. 



PACS numbers: 52.80.-s,51.50.+v,52.27.Aj,52.27.Cm 



Diffusion correction to the Raether-Meek criterion 



2 



1. Introduction 

1.1. Problem setting and review 

Emergence and propagation of streamer- like discharges are topics of current interest. 
Streamers play a role in creating the paths of sparks and lightning [1, 2] and 
in sprite discharges at high altitude above thunderclouds [3, 4, 5]. They are 
also used in various industrial applications [6], e.g., in corona reactors for water 
and gas treatment [7, 8, 9, 10], and as sources of excimer radiation for material 
processing [11, 12, 13], for a recent overview see [14]. 

In the present paper, we investigate the conditions under which a tiny ionization 
seed as, in particular, a single electron in a homogeneous electric held far from 
any electrodes grows out into a streamer with self-induced space charge effects and 
consecutive rapid growth. The critical length of time for this transition as a function 
of the electric held is usually described by the Raether-Meek criterion. We will 
confront current simulation results with the underlying assumptions of the Raether- 
Meek criterion, and then derive a diffusion correction to it. This correction can amount 
to a factor of 2 or more for transition time and length for certain parameters as we will 
elaborate below and summarize in Figs. 5 and 6. The consequences are particularly 
severe in non-attaching gases, where in low fields the diffusion can suppress streamer 
formation almost completely while the Raether-Meek criterion would predict streamer 
formation after a finite travel distance and time. An example of such an avalanche in 
extremely low helds is discussed in [15]. 

In many applications, discharges are enclosed by containers and electrodes; 
streamers then frequently emerge from point or rod electrodes, that create strong local 
fields in their neighborhood [16] and also influence the discharge by surface effects. 
On the other hand, in many natural discharges and, in particular, for sprites above 
thunderclouds [5], it is appropriate to assume that the electric field is homogeneous 
and metal electrodes absent. In this case, single electrons can create ionization 
avalanches that move into the electron drift direction. From those avalanches, single 
or double ended streamers can emerge, and we are interested in the prediction of this 
transition. For clarity, we call a spatial distribution of charged particles an avalanche, 
if the electric field generated by their space charges is negligible in comparison to the 
background external field; on the other hand, if the space charges of the distribution 
substantially contribute to the total field, we speak of a streamer. 

The critical field required for lightning generation is presently a topic of debate, 
in particular, whether thundercloud fields are sufficient for classical breakdown or 
whether relativistic particles from cosmic air showers are required [17, 18]. Different 
critical fields can be defined for different processes; e.g., in [16] a critical field for 
the propagation of positive streamer propagation is suggested that is valid after the 
streamers have emerged from a needle electrode. This held is certainly lower than the 
critical field for streamer emergence from an avalanche to be discussed here. 

Of course, dust particles or other nucleation centers can play an additional role 
in discharge generation in thunderclouds, but in the present paper we will focus on 
the effect of a homogeneous field in a homogeneous gas. This assumption corresponds 
to the classical experiments of Raether in the thirties of the last century [19]. 

Within the present introductory and motivating section, we first recall the 
common discharge model and present simulation results for avalanches and consecutive 
streamers that emerge from a single electron in a homogeneous held far from any 
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surfaces. Then we recall the Raether-Meek criterion; it suggests that the avalanche 
to streamer transition depends on the ionization rate a and gap length d through 
the dimensionless combination ad. We confront this criterion with our simulations 
and argue that the transition depends not only on the ionization coefficient times gap 
length but also on electron diffusion. Now numerical evaluations of the initial value 
problem for a large range of parameters, namely fields, gas types and densities, would 
be very tedious. However, we have succeeded in making analytical progress on the 
transition criterion. This has two major advantages: first, general expressions for 
arbitrary fields, gases and densities can be derived. Second, the result can be given 
in the form of a closed mathematical expression. These calculations and results form 
the body of the paper. 

1.2. Discharge model and simulation results 

In detail, we consider a continuous discharge model with attachment and local field- 
dependent impact ionization rate and space charge effects. It is defined through 
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where charged particles are present only in a bounded region, and the electric field 
far away from the ionized region is homogenous. Here n e , n + and n_ are the particle 
densities of electrons, positive and negative ions, and E and <& are the electric field and 
potential, respectively. The total field E is the sum of the background (Laplacian) field 
Ef, in the absence of space charges and the field generated by the charged particles E'. 
D e and v a are the electron diffusion and the electron attachment rate, respectively. 
The impact ionization coefficient on is a function of the electric field, as established in 
various books, and for our numerical calculations, we use the Townsend approximation 

a j (|E|) = a ex P (-fb/|E|) ) (5) 

in which a and E are parameters for the effective cross section. They depend on the 
ratio of background and normal gas density (N and No, respectively) as ckq oc (N/Nq) 
and Eq oc (N/N ) [20]. This scaling is equivalent to stating that the reduced electric 
field E/iV is the relevant physical variable for impact ionization processes. The positive 
and negative ions are considered to be immobile on the time scales investigated in this 
paper because avalanches and streamers evolve on the time scale of the electrons that 
are much more mobile due to their much lower mass. 

We consider the situation where a tiny ionization seed of the size of one or a few 
free electrons is placed in free space, i.e., within a gas far from walls, electrodes or 
other boundaries. If the externally applied field is sufficiently high, it will develop 
into an electron avalanche that will drift towards the anode. Eventually, the charged 
particle density in the avalanche will become so large that space charge effects set in 
and change the externally applied field. As a consequence, the interior of the formed 
very weak plasma will then be weakly screened from the external field while the field at 
the outer edges is enhanced. Depending on photo-ionization processes, then an anode- 
directed or a double ended streamer emerges from the avalanche. This evolution from 
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an electron avalanche to a streamer is illustrated in Fig. 1. Details on our simulations 
can be found in [21, 22, 23], here we only use them for illustration purposes. 

Fig. 1 shows essential features of the solutions that will be substantiated by 
quantitative analysis in the body of the paper. In the left column, an avalanche can 
be seen: the electron distribution (upper row) is Gaussian and spherically symmetric. 
The position of the Gaussian is determined by electron drift in the homogeneous 
background field, its width by electron diffusion. The ions (second row) are left behind 
(i.e., further down) and stretched along the temporal trace of the avalanche. The 
resulting space charge distribution (third row) is essentially a smooth dipole without 
much structure. Actually, these pictures are quite similar to the sketches of Raether. 
The electric field (fourth row) is essentially unchanged up to corrections below 1%. 
The current (lowest row) shows the same Gaussian structure as the electrons; it is 
dominated by electron drift /U e n e (R, t) in the homogeneous background field Ef, 
with a small diffusional correction. In the right column, a conducting filament is 
formed, and the streamer stage is reached. Electron and ion distribution show a 
similar long stretched shape. The space charges approach a layered structure, and the 
field ahead of the streamer is changed by these space charges by up to 40%. 

There is some freedom in defining the transition point from avalanche to streamer. 
In the body of the paper, we will argue that a maximal field enhancement of 3% ahead 
of the streamer, i.e., 

fc= max R |E(R t )|-lE & l fc = ^ 

l E b| 

is a decent measure for the transition. We will see that essentially up to this moment 
of time the total number of electrons in the avalanche grows exponentially in time, 
while in the streamer phase, the growth is slower. 



1.3. Review of critical field and Raether-Meek criterion 

Essentially two criteria have been given in the literature for the emergence of a 
streamer from a tiny ionization seed, one for the required background field and one 
for the required space and evolution time. The first one is a necessary lower bound 
for the background field: the electric field has to be higher than the threshold field 
.Efc where the impact ionization rate overcomes the attachment rate. The ionization 
level can only grow if the rightmost term in Eq. (1) is positive, hence if the effective 
ionization coefficient is positive, 

a(|E|) = ai(E|)-i/ o /(Ai e |E|)>0. (7) 
This determines the threshold field E k as 

He E k a.i{E k ) = v a . (8) 

The second criterion is known as the Raether-Meek criterion. It states that the 
total electron number must have reached the order of 10 8 to 10 9 for space charge 
effects to set in. If this number is reached by exponential multiplication of one initial 
electron within a constant field E&, this means that 

e a(|B|)d „ 1Q 8 tQ 1Q 9 } ( 9 ) 

where d is the avalanche length. In brief as a rule of thumb the criterion reads 

a(|E|) d w 18 to 21 according to Raether and Meek. (10) 
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Figure 1. The avalanche to streamer transition: numerical solution of the 
discharge model (l)-(4) for pure N2 in a uniform background field, under 
assumption of cylindrical symmetry. As N2 is a non-attaching gas, there is no 
formation of negative ions. The background electric field is directed in the negative 
z-direction and has a strength E b / {N / Nq) = 100 kV/cm where N is the actual 
particle density and No the particle density under normal conditions. Initially, a 
single electron is placed at a distance z=115 fim of the cathode (which in turn is 
placed at 2 = 0). 

Shown are the electron avalanche phase (left column), the transition to 
streamer (middle column) and the space charge dominated streamer phase (right 
column). The respective times are t= 0.225, 0.375 and 0.525 ps for N2 under 
normal conditions. 

From top to bottom: electron number density distribution; ion number 
density distribution; net charge density distribution (positive: thin lines, negative: 
thick lines) and equipotential lines (dashed); electric field strength (smaller 
than the background field: thin lines, larger: thick lines); current density 
j c = M,;En e + D e Vn e . 
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Let us first note that the same criterion has been suggested for quite different situations 
in the iiterature. In his original article, Meek [24] studies the emergence of a cathode 
directed (i.e., positive) streamer from an anode directed avalanche that has bridged 
a short gap. On the other hand, Bazelyan and Raizer [25] study the emergence of 
streamers in free space, i.e. far away from the electrodes. To estimate the field of 
the ions, Meek used the diffusion radius of the electron avalanche, and the ionization 
rate in the background field, however, the diffusion does not show up in his transition 
criterion. Bazelyan and Raizer on the other hand, neglect diffusion and base their 
criterion on the radius of the avalanche due to electrostatic repulsion. All authors 
[24, 25, 26, 27] assume the electron distribution to be spherically symmetric, on the 
other hand, they base their transition criterion on a total field screening, i.e., to k = 1 
in Eq. (6). In view of available simulation results like our Fig. 1, these assumptions 
are self-contradictory. 

Apart from these considerations on the history of the derivation, there are actually 
two major reasons to revise the Raether-Meek-criterion: 

1) The prediction that a parameter should be in the range of 18 to 21 (where authors 
seem to be willing to assume an even larger range of values to get consistency 
with experiment) is not very satisfactory and invites improvement. 

2) Diffusion has to be included into the model for physical as well as for mathematical 
reasons. Without diffusion, an initially concentrated electron package would not 
spread and it would create enormous fields within very short time as they are 
well known in the neighborhood of point sources. Indeed, diffusion decreases the 
electron density and the maximal fields while impact ionization increases it. In 
low fields, diffusion stays dominant for a long time and delays space charge effects 
and consecutive streamer emergence. It is therefore clear that the avalanche to 
streamer transition docs not only depend on multiplication rates, but also on 
the relative importance of diffusional spreading. This should provide a more 
quantitative transition criterion than the pure Raether-Meek criterion. 

1.4- Organization of the ■payer 

We will derive a diffusion correction to the Raether-Meek criterion through the 
following steps: in Section 2, the intrinsic scales of the problem with their explicit 
density dependence arc identified through dimensional analysis. In Section 3, we 
analyze the spatial distribution of the electrons during the avalanche phase and 
their contribution to the electric field; this gives a first approximate correction to 
the Raether-Meek-criterion. In Section 4, we approximate the spatial distribution 
of the ions and their contribution to the electric field. Electron and ion field are 
then combined to give the total change of the electric field during the avalanche 
phase. If this field becomes "substantial" (cf. Fig. 1 and Eq. (6)), we have found 
the avalanche-to-streamer transition. Finally, the analytical non-dimensional results 
are translated back to dimensional quantities, and we refer the reader interested in the 
final prediction only to Figs. 5 and 6 for the transition criterion in non-attaching gases 
and in air. These figures visualize the analytical criterion (51). Section 5 contains the 
conclusions, and an additional approximation for the ion generated electric field can 
be found in the appendix. 
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2. Dimensional analysis 

The Raethcr-Mcck criterion can be understood as a simple example of dimensional 
analysis. Dimensional analysis identifies general physical properties in terms of 
dimensionless numbers that are independent of a particular gas type or density. The 
physical importance of dimensionless numbers like the Reynolds number is well known 
in hydrodynamics, we here follow the same approach. 

In the light of dimensional analysis, the Raether-Meek criterion states that 
the effective cross-section a(|E|) has the dimension of inverse length, hence the 
dimensionless number a(|E|)d characterizes the gap length in multiples of the 
ionization length and therefore the exponential multiplication rate e ad . This number 
directly characterizes the total number of electrons in an avalanche that started from 
a single free electron. However, this is not the only dimensionless number in the 
problem, a second one is the dimensionless diffusion constant 

that plays a distinctive role in the avalanche to streamer transition as well as it 
determines the width of the electron cloud. Note that this dimensionless diffusion 
constant is related to the electron temperature as D e /fi e — ksT e where ks is the 
Boltzmann constant. The electron temperature T e actually can be defined through 
this relation, even if the electron energy distribution is not Maxwcllian in the presence 
of strong electric fields. Furthermore, D depends on ao/Eo, the two parameters 
characterizing the Townsend coefficient, given by Eq. (5), for the specific gas and 
density. 

For the general performance of dimensional analysis, we refer to earlier articles 
[28, 29] and here only state the results: lengths are measured in units of 1/ao, fields 
in units of Eq, velocities in units of vo = fi e Eo and time consistently in units of 
to = 1/ (ao/i e Eo) — hence diffusion should be measured in units of /i e Eo/ao as done 
in (11). The natural scale for the particle densities follows from the Poisson equation, 
n a = e a a a E a /e. 

The parameters a>o, He, D e and Eq depend on the ratio of the background gas 
density N and the gas density under normal conditions N . Using parameters as 
in [20, 30, 31, 32], the characteristic scales are for N 2 : 

, 1 kV N 

^ = 380 ^-)-, D e = 1800 ^-)-, 

Vs N/No 1 s N/No 1 (i2) 

„ ,„ 7 cm N y ' 

vq = 7-W — , t = 3ps— , 

s , N ° 

no = 4.8 • 10 14 ^ 

cm 3 V^o, 

and the dimensionless diffusion constant is D = 0.1. Notice that the characteristic 
velocity scale is independent of pressure, in agreement with measurements of streamer 
velocities at different pressures. Notice furthermore, that it directly follows from this 
analysis that the relevant physical parameter is the reduced electric field E/N. 
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Dimcnsionless parameters and fields are introduced as 

R t 

r=-r, t=— , i/ = v a t , 

l to 

n e n+-n_ E 

n n Eo 

which brings the system of equations (l)-(4) into the dimensionless form 

d T <7 = D\7 2 a + V(£<r) +f(\£\,v)a, (14) 
d rP = f(\£\,v)a, (15) 
-V 2 (f> = p-a £ = -V<t>, (16) 

where the operator V is taken with respect to r and where f{\£\, v) is the dimensionless 
effective ionization rate, 

fm^) = ^ m =\S\e-^-,. (17) 

He -fc/O CIO 

It is remarkable that the density of positive and negative ions n± enters the equations 
only in the form of the single dimensionless field p oc n + — n_ . This is clear in the 
case of the Poisson equation, but holds also for the generation term proportional to 
f(\£\,v). This coefficient accounts for the production of free electrons through impact 
ionization and for the loss of free electrons due to attachment. 

We neglect the effect of photoionization as its rates are typically much lower 
than impact ionization rates; it does not contribute significantly to the build-up of 
a compact ionized cloud where eventually space charge effects will set in (quite in 
contrast to its distinct role in positive streamer propagation). 



3. Electron distribution and field 



We derive the transition as follows: We assume that an avalanche starts from a single 
electron and follows a transition as shown in Fig. 1. In the calculation we neglect 
space charge effects on the evolution of densities, but we do calculate the additional 
electric field generated by the space charges. If this field reaches a relative value of 
k=0.03 — this value will be motivated in Section 4.3 — , space charge effects are not 
negligible anymore, and the transition to the streamer is found. 

The electric field generated by space charges has one contribution from the 
electrons a and another one from the positive and negative ions p. In the present 
section, we calculate the field of the electrons, in the next section, we will include the 
field of the ions. 



3.1. The electron distribution: a Gaussian 

We write the single electron that generates the avalanche as a localized initial density 

cr(r,T = 0) =p(r,T = 0) =cr <5(r-r o ) (18) 

and consider its evolution under influence of a uniform field £b = —£b& z , where e z is 
the unit vector in the z direction and £b = \£b\ is constant. A single electron is written 
as a ((-function n e (R) oc 5 3 (R — R ) in physical units where the spatial integral over 
the electron number density 

N e {r) = /d 3 Rn e (R), (19) 
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of course, should be unity N e (0) = 1. According to the last section, this corresponds 
in dimensionless units to cr = l/( n o^o) which is 1.7 • 10~ 4 N/N for nitrogen. We 
will use cr = 10~ 4 in the sequel. We emphasize, however, that the theory will be 
developed for an arbitrary value of ao- 

During the avalanche phase the electric field remains unaffected by space charges, 
so that the continuity equations for the charged particles (14)-(15) can be linearized 
around the background field, 

d T <j = DV 2 a + £ b - Va + af, (20) 

T p=af, (21) 

where / = f(£ b , v). 

For the initial condition (18), the electron evolution according to Eq. (20) can be 
given explicitly as [20] 

tlcr s e -(r-v a +£ b T) 2 /(iDr) 

«M=a ef^ _____ . (22) 

it has the form of a Gaussian package that drifts with velocity —£b, widens diffusively 
with half width proportional to ^/ADt, and carries a total number of electrons 
aoe^ Sb ' v ' T . (If the initial ionization seed consists of several electrons in some close 
neighborhood, the Gaussian shape is approached nevertheless for large times due to 
the central limit theorem.) 

Integrating Eq. (22) over the entire space shows that the total number of electrons 
grows as N e (r) — a n ile^ T (if we start with a single electron). On the other hand, the 
maximum of the electron density is reached at the center of the Gaussian r = r — £ b r 
and has the value 

<7max(T) = maxa(r,r) = ^ t)3/2 , (23) 

hence it first decreases until r = 3/(2/) due to diffusion and then increases due to 
electron multiplication. At this moment of evolution, generation overcomes diffusion. 

The axial electron density distribution for a background field of £ b = 0.25 at 
r = 2000 (for N 2 this corresponds to a reduced electric field E b (N /N) = 50 kV/cm 
and t—6 ns) is illustrated in the upper panel of Fig. 2. The analytical solution (22) of 
the linearized continuity equation (20) is compared to a numerical evaluation of the 
full nonlinear problem (14)-(16). The excellent correspondence between the solution 
of both the linearized and the nonlinear problem shows that, at this time, space charge 
effects are negligible, so that the electrons still are in the avalanche phase. 



3. 2. Exact result for the electron generated field £ a 

While density and field of the ions can only be calculated approximately and will be 
treated in the next section, the electric field £ a generated by the Gaussian electron 
package can be calculated exactly. 

The main point is that the electron density distribution (22) is spherically 
symmetric about the point r — £ b T. The electric field £ a (s 7 r) at the point 

s = r - r + £ b r (24) 

can therefore be written as £ a {s, t) = — £ a {s, r)e s , where e s is the unit vector in the 
radial s direction. Its magnitude can be computed with Gauss' law of electrostatics 
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Figure 2. Analytical results (solid lines) for the linearized problem compared 
to a numerical evaluation (symbols) of the full nonlinear model (14)-(16) in a 
background field £ b = —0.25 e z , and D = 0.1, v = and crn = 10~ 4 . The 
time of the snapshot is r = 2000. The initial condition is located at rn = 50 e z . 
Upper panel: the electron density distribution (x) and ion + density distributions 
along the axis of symmetry r=0, as computed by a numerical simulation of the 
nonlinear model; the solid line is the analytical solution (22) of the linearized 
problem. Lower panel: the axial component of the numerically computed space 
charge field (i.e. generated by the electrons and the ions) £' = £„ + £p = £ — £b 
(thin line) and the analytical result for the field £ CT generated by the electrons 
only (thick line). The field induced by the ions is discussed in the next section and 
not accounted for in £ a . We remark that this situation corresponds qualitatively 
to the leftmost column of Fig. 1. 



(in the same way as the gravitational force field of a spherically symmetric mass 
distribution). It uses the fact that the field at radius s is determined by the total 
charge inside the sphere of radius s, and independent of charges outside this radius, 
as long as the distribution is spherically symmetric. It yields 

with 

F{x) = — — / y 2 e y dy = — , 26 

x Jo x 2 0r x 

where erf is the error function. 
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The spatial maximum of the field strength £ a is determined by the maximum of 
F(x); evaluating F'(x) — shows that it is located at an x such that 

-?=(x + x 3 )e- x2 =erf x. (27) 
v 7r 

Solving this equation numerically leads to a position of the maximum of about x ~ 1 
(which is the radius at which the Gaussian electron distribution has dropped to 1/e 
of its maximal value) and to the value F(l) ~ 0.4276. The spatial maximum of the 
electron generated electric field strength becomes 

^ (T) "il^ F(1) ' (28) 

it is located on the sphere parameterized through 

|r-r -5 6 r| ~ Vlth 1 . (29) 

In the original cylindrically symmetric coordinate system (r, z) , the axial field 
component is directed in the negative z-direction, i.e. in the same direction as the 
background field, "ahead" of the electron cloud (z > zo + £br) as is illustrated by the 
solid line in the lower panel of Fig. 2. Combining this with Eq. (29), we find that the 
maximal field strength \£ b + S a \ and its location are 

mzx\£ b + £ a \ = \£ b + £ a \(r m ,T)=£ b + £™ ax (r), (30) 



r m (r) ~ (z + £ b r + VlD^)e z . (31) 



3. 3. A lower bound for the transition 

Since the avalanche to streamer transition takes place when space charge effects start 
to affect the electric field, we choose to base the criterion for the transition on the 
maximal relative field enhancement fc(r) defined in Eq. (6), which for the dimensionless 
field simply reads 

= max r |£(r,T)M^ (32) 
\£b\ 

Here £ = £ b + £ a + £ p is the total electric field, £ a and £ p being the fields of the 
electrons and the ions, respectively. We will show in the next section that kt = 0.03 
is an appropriate estimate for the maximal relative field enhancement at the mid gap 
avalanche to streamer transition. At lower values of k, space charge effects can be 
neglected, whereas at higher values the dynamics of the electrons are nonlinear and 
the full streamer equations (22)-(16) have to be solved. 

As a first estimate for the space charge field, and thereby for the avalanche to 
streamer transition, we compute the field generated by the electrons only and neglect 
the ion field. This is a decent approximation, as the lower panel in Fig. 2 shows. 
Actually, the magnitude of the monopole field £ a ahead of the electron cloud is 
an upper bound for the magnitude of the field created by the dipolc of electrons 
on the one hand and the positive charges left behind by the electron cloud on the 
other hand. Therefore, the maximal relative field enhancement due to the electrons, 
fco-(r) = £™ ax ( T )/£b, exceeds the transition value after a shorter travel time t ct and 
distance then the genuine relative field enhancement k(r) of Eq. (32). Hence, T a is a 
lower bound for the time r a ^ s of the avalanche-to-streamer transition. 
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The lower bound r CT for the transition can be expressed through Eq. (28) as 

/tct - ]n(£ b T a ) ~ In ~^y~ - (33) 

As travel time and travel distance are related through the drift velocity £(,, ^)t ct 
is found to be identical to (a(|Et,|) — v a j p e Eif)d a in dimensional units where d a = 
\x e E\,t a is the avalanche travel distance. In dimensional quantities, Eq. (33) takes the 
form 

i Tp i\ v a \ 167rfcl0 4 D e a N 

For a non-attaching gas [y a = 0) at atmospheric pressure under normal conditions 
with dimensionlcss diffusion comparable to nitrogen, inserting the numerical values 
for the parameters, we obtain 

a(|£ 6 |)d CT -ln(a ck)~ 9.43. (35) 

/ being a growing function of \£b\, Eq. (33) shows that the larger the field, the earlier 
the transition takes place, which is in accordance with Meek's criterion. On the other 
hand, the second term on the right hand side of Eq. (34) depends on the diffusion 
coefficient in such a way that diffusion delays the transition to streamer, as expected. 

The solution a(\Eb\)d a for N 2 at atmospheric pressure is shown in the dash-dotted 
line of Fig. 3, where it is compared to a numerical evaluation of the transition time 
(circles). The latter have been obtained through a full simulation of the continuity 
equations (14)-(15) together with the Poisson equation (16) [29, 21] that was also 
used to generate Fig. 1. Though the qualitative features of the transition time are 
well reproduced, this figure shows that the underestimation of the transition time is 
significant, and that it is necessary to include the field of the ion trail left behind by 
the electrons. 

4. Ion distribution and field 

4-1. Exact results on the spatial moments of the distributions 

To get a more accurate estimate for the avalanche-to-streamer transition, the field 
generated by the positive and negative ions has to be included. In the case of the 
ion distribution, closed analytical results cannot be found, in contrast to the electron 
distribution (22). However, arbitrary spatial moments of the distribution 

f O o d 3 r 

(°)p = —r n . wherc ° = zn or ( 36 ) 
J p d 6 v 

can be derived analytically. Here z is the direction of the homogeneous field £ b and 
r is the radial direction. First, the evolution equation (15) for the ion density p is 
integrated in time and the analytical form (22) for <j(r, r) is inserted. As / — f(\£b\, v) 
is constant in space and time one finds 

p(r,r)-p(r,0)= (37) 

-(z-z -£„T') 2 /(4Dr') p -r 2 /(ADr') 



fo- f dr' ef T ' 
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Figure 3. The transition time fr (equivalent to the travel distance ad) as a 
function of the background electric field for ao = 10~ 4 , v = and different values 
of D. Solid line: computed with Eq. (51) for D = 0.1 (thin line), 0.3 (medium thin 
line) and 1 (thickest line); dash-dotted line: computed with Eq. (33) for D = 0.1; 
symbols: numerical evaluation for D = 0.1. Obviously, the approximation (51) 
fits the full numerical results very well. 



Here the initial perturbation is located at z on the axis r = 0. The moments (36) can 
now be derived from (37) by exchanging the order of spatial and temporal integration. 
In particular, one finds 

J pd 3 r = o~ o e /T , (38) 

J z p d 3 r = ao e }T + £ b r - —JJg^j > 

and higher moments can be calculated in the same way. For the moments of p in the 
axial direction, this gives 

(z) P =z +£ b (r--\ +0(e~f T ), (39) 



zo+£ b (r-j)+0(e-t T ), 

<A> (j) 2 + 2D ( T -j) + °( e ~ fT ) 



(40) 



The second moment of p in the radial direction is 

{r 2 ) p = 2 D (r--^+0(e-^). (41) 

For comparison, the moments of the Gaussian electron distribution (22) are easily 
found to be 

{z)a =z + £ b T, (42) 

(z 2 ) a - (z)l - 2Dt, (43) 

(r 2 ) a = 2Dt. (44) 
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4-2. Discussion of the moments 

Let us now interprete these moments. A first moment of a spatial distribution gives 
its center of mass. For the second moment, the cumulant 

{**)%•■=({* -{*)*?) x = {*%-{z)l x = a,p. (45) 

measures the quadratic extension from the center of mass. As the center of mass lies 
on the axis, for the radial extension the distinction between second moment and its 
cumulant need not be made. 

The moments for the electrons (42)-(44) have a simple structure: the center of 
mass of the electron package is located at z = z^ + E^t, and the package has a diffusive 
width \/1Dt around it, both in the forward z direction and in the radial r direction. 

The ion cloud shows a more complex behavior; it is evaluated close to the 
avalanchc-to-streamer transition where fr = ad = O(10), therefore the terms of 
order e~^ T are neglected. 

First it is remarkable that the center of mass of the ion cloud (39) shifts with 
precisely the same velocity as the electron cloud though the ion motion is neglected 
while the electrons drift, therefore the ion center of mass is at an approximately 
constant distance £&// behind the electron center of mass. This distance 

} (Lb) a(Eb) 

corresponds to the dimensional ionization length \/a(Eb). 

The quadratic radial width of the ion cloud 2D(t~ 1//) is 21?// smaller than the 
one of the electron cloud. This is related to the fact that the electron cloud also was 
more narrow at the earlier times when it left the ions behind. The ion cloud is more 
extended in the z direction. More precisely, its length is £ a larger than its width. This 
comes from the ions being immobile, therefore a trace of ions is left behind by the 
electron cloud. Moreover, it can be remarked that the difference between quadratic 
width and length of the ion cloud is given by the same ionization length £ a as the 
distance between the centers of mass of the ion and the electron cloud. We refer to 
the left column of Fig. 1 for the illustration of these density distributions. 

4-3. An estimate for the transition 

One can assume as in [25] that the ions have a spatial distribution similar to the 
electrons, thus a Gaussian with the same width as the electron cloud, but centered 
around (r = 0, z = (z) p ): 

e -[(z-{z) p ) 2 +r 2 ]/(4Dr) 

Pl (r,z,r)=a e^ { ^ Dt)V 2 ■ ( 47 ) 

In this approximation, the total electric field becomes: 

oft 

£i(r,z,r) = £ b - 



F I 1 — + F S ' 8 



KVaD^J I s^l \V4~DtJ K 



(48) 



16tt Dt 
where 

s x = r - (z) x e z for x = p, a (49) 
are the distances to the electron and ion centers of mass. 
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The maximum of the field £\ can not be computed analytically. However, in 
Fig. 1 and in the lower panel of Fig. 2, it can be seen that this maximum is located 
on the axis ahead of the electron cloud, and that the location of the maximum of 
the total field and that of the electron field nearly coincide. This can easily be 
explained physically: the total field is the sum of the fields induced by the electrons 
and by the ions. Its maximum is located just ahead of the electron cloud, where the 
electron field varies rapidly, while the field contribution of the ions is smoother and 
smaller as we are interested in its contribution further away from the center of the ion 
distribution. Therefore the maximum position of the total field is essentially identical 
to the maximum position of the electron field. This justifies our approximation to 
evaluate the field E\ at the maximum position r m of \£ b + S a \ as defined in Eq. (31). 
The maximum of the electric field can thus be approximated as: 



£™ ax (r) ~ &(r = 0, z = z + £ b r + VId7, r) 

ofr ( I 02 

F(l)-F[l + 



£ b + a ° e 



IQttDt 



ADt 



(50) 



Then £™ ax — £ b = k£ b implies for the transition time n : 

fn - H£ b n) - In — , = m (51) 

where F(x) is defined in Eq. (26). The argument of the logarithm in the third term on 
the left hand side is larger than 1, therefore this criterion gives a later transition time 
than that based on the field of the electrons only. This is what we expect considering 
that the ions tend to reduce the field of the electrons, thus the effect of space charge. 
The correction given by the ion field is a function of the ratio of the ionization length 
£ a and the diffusion length y/2Dr. At early times, this ratio goes to infinity, and the 
correction given by the ion cloud is negligible. However, at later times, the correction 
becomes more significant. 



4-4- The analytically approximated transition criterion compared with numerical 
results 

We now compare again our analytical results for the linearized problem to the outcome 
of numerical simulations of the full nonlinear model (14)-(16). 

In the upper panel of Fig 4 the evolution of the maximal electron density as a 
function of /(|£f,|)r is shown. Numerical and analytical solutions coincide during the 
avalanche phase, but deviate eventually. This enables us to estimate the moment at 
which the space charge effects set in, and thus when the streamer regime is reached. 
In the lower panel of Fig 4 the evolution of the maximal relative field enhancement is 
considered. Looking at the simulation results (the solid lines), we see that k = 0.03 
gives a good estimate of the transition time. 

The approximation (50) for the maximal field has now become much better than 
the previous approximation (28) based on the electron cloud only. Indeed, for e.g the 
case of £ b = 0.5 (corresponding to the medium thick lines), the numerically computed 
field (solid line) reaches the transition value ((£ nU m — £ b ) = 0.03£& at fr w 14. 
When only the field of the electrons is taken into account, this value would already be 
reached at fr w 12.6, while the correction based on the approximation of the ion cloud 
leads to a transition time of fr w 13.9. The correction becomes especially important 
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Figure 4. Comparison of the analytical approximations and simulation results 
of the full non linear streamer equations (14)-(16), for various background electric 
field strengths £{, ; thin line: £\y = 1, medium thick line £ b = 0.5 and thick line 
£ b = 0.25. 

Upper panel: the evolution of the maximal electron density as a function of /(£;,t) 
as computed within the full nonlinear 2D model (solid lines) and as given by the 
analytical solution (23) of the linearized problem (dash dotted lines). 
Lower panel: the evolution of the maximal electric field enhancement k = (£ m ax — 
£b)/£b as a function of f(£ b r). Solid lines: numerical solution of the full nonlinear 
2D model; dashed dotted lines: only the field of the electrons is accounted for in 
the analysis, see Eq. (28); dashed lines: analytical approximation (50) of the total 
field. 



at higher fields. In low fields, the approximation of the ions shows somewhat larger 
deviations. We notice that the analytical approximation p\ is narrower and higher 
than the genuine one, and therefore leads to an overestimation of the field generated 
by the ions inside the ion cloud. For an even more accurate estimate of the total field 
between the electron and the ion cloud we refer to Appendix 2, where it is also shown 
this will not lead to a significant improve in the estimate of the maximal field ahead 
of the electron cloud. 

In Fig. 3 we compare the transition times given by Eqs. (33) and (51) with 
numerically evaluated transition times. It shows that the approximation of similar 
electron and ion distributions leads to a very good approximation of the transition 
time. This figure also illustrates that the transition time fr depends strongly on the 
electric field, and increases for smaller fields. Moreover, looking at the transition time 
for higher diffusion coefficients, it is seen that diffusion tends to delay the transition 
to the streamer regime. This can be expected, since diffusion will tend to broaden the 
electron cloud, thereby suppressing space charge effects. Depending on the external 
parameters, the value of ad at the time the transition takes place can vary by a factor 
two or more. 
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4-5. The final results on the transition criterion 

The transition time approximated by Eq. (51) as a function of both background 
electric field and diffusion coefficient is visualized as the 3-dimensional in Fig. 5. This 
figure shows that the Raether-Meek transition criterion, that stated that fr takes 
an approximately constant value of 18 to 21, corresponds to the case of relatively 
high diffusion and background field. However, realistic values of D are smaller than 
unity, and a background electric field higher than 2 also leads to unrealistic values. 
So in the parameter range of real experiments, the correction given by Eq. (3) on the 
Raether-Meek criterion can not be neglected. 




deviate from the Raether-Meek criterion (10). 

We no discuss the particular example of an electron avalanche in (dry) air, 
for which different coefficients have to be used than in N2. In particular, the 
ionization length and field in air are given by [20] ao = 0.87/Ltm(iV/iVo) and Eq = 
277.4kV/cm • N/Nq- For the values of the mobility and the diffusion coefficient of 
the electrons as a function of the electric field we use experimental values as well as 
numerical solutions of the Boltzmann equation (see Appendix A). Inserting those in 
Eq. (51), we can compute the value of a(|E|)d at the transition for different background 
fields, showed in Fig. 6. At large fields, the value of ad at transition saturates towards 
16, and grows asymptotically as the reduced field approaches a value of ±27.7 kV/cm. 
At even lower fields attachment overcomes electron impact ionization, and a single 
electron can not generate a streamer. Large values of ad as in Fig 5 are not found, as 
electron attachment limits the emergence of streamers in low fields (see Eq. (8)). So 
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Figure 6. The value of ad when the electron avalanche has traveled a distance 
d and has reached the point at which the transition to streamer takes place. The 
x-axis shows the dimensionless electric field. The line at £ ~ 0.01 indicates the 
field below which attachment overcomes ionization, and hence there is no electron 
multiplication, o: experimental value for a(E) [33] and D/mu e [34]; □: values 
from the Boltzmann solver [5] . 



for air, ad drops from 21 to 15 with growing field. 
5. Summary and conclusions 

Recent simulations (see Fig. 1) have shown that an electron avalanche turns into 
a streamer when the field enhancement due to space charges is about 3%. In this 
paper, the theory behind the commonly used Raether-Meek criterion is reviewed, as it 
assumes a linear behavior of the electrons (i.e. space charge effects can be neglected), 
which is in contradiction with the requirement of the space charge field to be in the 
same order of magnitude as the background field for the transition to occur. 

A dimensional analysis has been carried, enabling us to identify the characteristic 
length scales, which are a function of the neutral gas density. In particular, the 
dimensionless quantities a(£ b )d and D = D e a /(^ e E ) have been extracted from the 
problem. The first gives the distance d in multiples of the effective ionization length in 
the background field, while the latter gives the ratio between diffusive and advective 
transport of electrons. The continuity equations for the positive and negative ions (2)- 
(3) reduce to one single equation (15) holding both positive and negative ions after 
rescaling, making the further analysis valid for both attaching gases like air or non- 
attaching gases like N 2 or Ar. 

The avalanche regime was identified as the phase during which space charge effects 
are negligible. This implies that the problem can be linearized around the background 
field, making it well-suited for analytical treatment. Indeed, for an electron avalanche 
evolving in a homogeneous background electric field, a closed analytical expressions 
exist for the density distribution of the electrons. Comparing this analytical solution 
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of the linearized problem to the results of a numerical simulation of the full nonlinear 
problem, it could be concluded that the transition to streamer takes when the maximal 
relative field enhancement k has reached a value of approximately 3%. 

We have shown that the electric field of the electron cloud during the avalanche 
regime can also be described by a closed expression. This led to the derivation of 
a lower bound for the avalanche to streamer transition (33). The estimate of the 
transition time has been improved by taking into account the field of the ions for which, 
in contrast to the electrons, no closed expression exists. However, the contribution of 
the ions to the maximal relative field enhancement can be well approximated, leading 
to an analytical estimate of the avalanche to streamer transition (51). 

The transition distance ad strongly depends on diffusion D and on the background 
electric field. For high fields, the transition time saturates towards ad ~ 15. On the 
other hand, for low fields, when the processes are diffusion dominated, the avalanche 
lasts longer. We remark that the striations observed in [15] are generic for atomic 
gasses with essentially only elastic and ionizing collisions, i.e. with very few inelastic 
processes [35]. 

In air, attachment limits the emergence of a streamer in low fields (see Eq. (8)). 
In this case, ad at transition is in the range of 16 (for high background fields) to 21 (for 
fields approaching Eh). It is remarkable that in the end, due to attachment cut-offs, 
Meek's criterion performs quite good for the emergence of streamers in free space. 
In non-attaching gases like N2 or Ar, the correction on Meek's criterion, that stated 
that ad ~ 18 — 21, becomes important at low fields. There the relatively strong 
diffusion delays the transition to streamer considerably. We emphasize that the use of 
dimcnsionless quantities enable us to translate the criterion given in (51) to any given 
neutral gas type and density. Evaluating the characteristic scales for these conditions, 
dimcnsionless field and diffusion can be computed, and the value of ad at transition 
can be computed from Eq. (51) or read from Fig. 5. Actually, Fig. 5 can also be used 
for attaching gases, as long as the ionization threshold field Ek is accounted for. 

The analytical models presented in this paper give a useful tool to describe 
the streamer formation. We stress that our criterion for the transition is based 
on a significant contribution of the space charges on the background electric field. 
Our analysis fully relied on the linearization of the streamer equations around the 
background field. The nonlinear streamer propagation is the subject of other studies. 
In that phase the space charges and electric field strongly interact, and the analytical 
study of such streamers [36] is far more difficult than the analysis of the linear 
avalanche phase. 
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Appendix A. Mobility and diffusion coefficients of electrons in air 

To compute the transition time in air, we use values of the electron mobility and 
diffusion coefficient found in literature. In the left plot of Fig. Al measured and 
calculated values of are given, as well as the fit = a exp(— E /\E\) 

are shown. The experimental values have been found in the survey of electron 
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swarm data by Dutton [33]. The computed values are the solution of the Boltzmann 
equation and have been taken from [5]. Also, the empirical approximation of the 
ionization coefficient as a function of the background field as given by [20] is shown, 
a(|E|) = a exp(-E /\E\) with a = 0.87pm(N/N ) and E = 277.4kV/cm • N/N . 

The values of D e j ' p e as a function of the reduced electric field are given in the 
right plot of Fig. Al. Again, computed values from [5] are shown, as well as measured 
values found in [34] . The value of the diffusion coefficient as a function of the electric 
field can obviously easily be derived from these figures. 




Figure Al. The ionization coefficient (left) and ration of electron diffusion 
and mobility (right) in air, as a function of the reduced electric field. o: 
experimental measurements (the values for a are taken from [33], the values for 
D/fi e from [34]. □: solution of the Boltzmann equation [5]. The solid line shows 
the ionization coefficient following the empirical formula (5) given in [20], with 
a = Q.87fj,m(N/N Q ) and E = 277.4kV/cm ■ N/N . 



Appendix B. A more accurate approximation for the ion density 
distribution 

The approximation for the ion distribution pi in Sect. 4 leads to a relatively good 
approximation for the transition time in the case of a mid-gap transition. However, the 
real spatial distribution of ions is more narrow in the r-direction, and can be wider and 
asymmetrical in the z-direction. In this appendix we present another approximation 
for the ion distribution, which will lead to a better overall approximation of the electric 
field, and of the self field induced by the ion trail. The price however to pay for this 
is a much more complicated analytical expression for the density and the field. 

A better approximation for p would then be an ellipsoidal Gaussian distribution 
centered around (r = 0,z = (z) p ) with width (z 2 ) c p = (z 2 ) p — (z) 2 and (r 2 ) c p = (r 2 ) p in 
the z- and r-direction, respectively. The height of this Gaussian should be such that 
the total amount of ions at time t is still equal to aoe?*. The appropriate expression 
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Figure Bl. The ion density (upper figure), total charge density (middle figure) 
and electric field (lower figure) on the axis, computed with Eq = 0.25, at r = 2000. 
The solid lines give the numerical solution, the dash-dotted lines the solution 
corresponding to p\ and the dotted lines to p2 . 
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for the ion distribution is: 

P(r,z,t) = ^g^ e- 2 /( 2 ^)-(^-^ P )V(^) (B.l) 

However, as far as we know, no closed analytical expression is known for the field 
of such an ellipsoidal Gaussian charge distribution. So instead, we take a spherical 
Gaussian distribution with the same height as the one defined in Eq. (B.l): 

P2(r, z, r) = e-^-^fVVSl) > (B . 2) 



where 



p 



= ^2D(r-j)^2D(r-j)+llj . (B.3) 
The electric field induced by this ion distribution is: 

where s p is defined in Eq. (49). 

In Fig. Bl we compare the densities and fields given by the numerical solution 
and pi and pi. It shows clearly that the approximation p 2 does not give a better 
approximation of the field ahead of the electron cloud. This can be explained by the 
fact that, the region ahead of the electron cloud does not contain any ions, so that the 
field induced by the ions is only a function of the total number of ions, which is the 
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same in both p x and p 2 - On the other hand, inside the ion cloud the approximation 
is much better. Therefore, evaluating the electron and ion densities with Eqs. (22) 
and (B.2) and their fields with Eqs. (25) and (B.4), at the transition time Ti given by 
Eq. (51), will give a good approximation of the status of the process at the time that 
streamer regime is entered. 
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